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I. INTRODUCTION. 



The evolution of the matter formed in ultrarelativistic heavy ion collisions is usually 
described within hydrodynamic and transport models. The former allow one to incorpo- 
rate the complicated evolution of the hot and dense system (fireball) at the possible phase 
transitions encoded in corresponding equation of state, while the latter make it possible 
to evaluate particle spectra taking into account the nonequilibrated character of their for- 
mation even at the early evolution stage. However, in realistic conditions, neither of these 
approaches nor their hybrids are able to reproduce particle spectra and correlations in A+A 
collisions. At present there is no scenario of the evolution allowing one to describe all the 
salient features of hadronic observables, including the interferometry or correlation radii re- 
lated to the space-time scales of the fireball. Particularly, the observed approximate energy 
independence of the radii is considered as the interferometry puzzle. 

It was noted Jl| , based on an analysis of the exact solution of the nonrelativistic Boltz- 
mann equation [2j, y], that under specific conditions which could be realized at the moment 
when fireball starts to decay, e.g., at the end of the pure hydrodynamic stage of the mat- 
ter evolution, spectra and interferometry radii may become "frozen" at this early moment 
despite the particle emission and collisions between hadrons last for a fairly long period of 
the system decay. It was also proposed [l( that this peculiar theoretical result can be used 
to test numerical transport models of A+A collisions and this is one of the aims of this 
paper. Another aim is a study, within kinetic approach, to what extent is the conservation 
of spectra and interferometry radii violated during the system evolution in the case of more 
realistic initial conditions. This study is particularly interesting in the context of recently 
observed [4] approximate conservation of interferometry volume during isentropic and chem- 
ically frozen hydrodynamic evolution of hadron-resonance gas and, it can be quite useful to 
solve the interferometry puzzle. 

The analysis of above phenomena has been done within Universal Kinetic Model (UKM) 
0|, a Monte Carlo event generator, realized as a numerical code written in object oriented 
C++ language. The UKM gives a possibility to model kinetic processes inside a particle 
system applying a suitable numerical algorithm. This model is a universal one, since there is 
a possibility to choose and add particle system (e.g., system of hadrons, system of partons, 
etc.) and to choose and add a particular numerical algorithm, i.e., to choose and add different 
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rinetic modeler. In our studies, we have used a cascade version of the UKM code, UKM-R 
I written under the ROOT Sy ste m Q. We ase the sa m e — eascade a lg o ri th,„ a8 
implemented in the Quark Gluon String model 

The paper is now organized as follows. In section |H] we relate particle spectra and two- 
particle correlation function with distribution and emission functions. Section ITTT1 deals with 
the theoretical aspects of the spectra formation, including a possible duality between kinetic 
and hydrodynamic approaches. Section |IV| is devoted to a brief UKM-R description. In 
section IS we discuss the known exact solution of the nonrelativistic Boltzmann equation 
for an expanding fireball with the spherically symmetric Gaussian initial conditions. This 
solution is used in section IV II to test the UKM-R. Section IV 111 deals with the influence of 
anisotropic initial conditions and relativistic effects on the final spectra and interferometry 
radii. The results are summarized in section IVIIII 



II. DISTRIBUTION AND EMISSION FUNCTIONS. 

The inclusive particle spectra are related to the averages of the products of creation and 
annihilation operators, a+ and a p , calculated at sufficiently large time t out or, generally, on an 
asymptotic hypersurface cr out : t = t out (x), guaranteeing the absence of further interactions. 
Thus 

d 3 N d 6 N 
P °~dFp ^ U ^ = W ^' P ° lP ° 2 rf3 pi #p 2 = nfa'P*) = ( a pi a t a Pi a p*)- i 1 ) 

In the case of independent (chaotic) particle production and subsequent system kinetic evo- 
lution, one can express the average of the product of four operators through the irreducible 
averages (a* ay), similar to Wick theorem for thermal systems. For identical spin-0 bosons 

n(pi,p 2 ) = «Si)( a P2 a P2) + («p 1 « P 2)( a p 2 a Pi)- ( 2 ) 

In kinetic approach, the irreducible averages are expressed through the integrals or Fourier 
transforms of the distribution function f(x,p) on the hypersurface cr out : 

(a+a p ) = / d 3 a IM (x)p f "f(x,p), «a P2 ) = / d^a^p^e^ f{x,p), (3) 

^°"out "'"but 

where, in the nondiagonal average in the second integral, p = (pi +p 2 )/2 is an off- mass-shell 
four- momentum and q = pi — p 2 . The actual choice of £ out does not matter due to the 



current conservation for free streaming particles and, due to the fact that, as a consequence 
of the relation qp = 0, the four-vector product qx is independent of the position x on the 
trajectory of a free off-mass-shell particle with four-momentum p (i.e., qx is invariant with 
respect to transformations t — > t + St and x — ► x + (p/p°)St). 

The distribution function f(x,p) is closely related with the current emission function 
S a (x,p) on a hypersurface a : t — i(x) since it collects contributions of the emission or 
scattering points x = {t, x — (p/p°) (t — t)} starting from those free particles with the three- 
velocity p/p° reach a point x = {t, x}. Particularly, choosing the hypersurface a : t — const, 
one has 

p°f(x,p) = J d 4 xS 3 (x - x + (p/p°)(t - t)) S a (x,p). (4) 

Inserting Eq. (jlj) into Eq. © and using the equality qx = qx, one can rewrite the asymptotic 
operator averages as four- volume integrals or Fourier transforms of the asymptotic emission 
function S(x,p) = S a ° ut (x,p) 0|: 

(a^ap) = J d A xS(x,p), (a^Opa) = J d 4 xe~ iqx S{x,p). (5) 

Using Eqs. (HJ-©, one can write the two-boson correlation function (CF) as 

n{ Vl ,p 2 ) = = 1 + \HP,q))? = 1 + (cos(gx 12 ))', (6) 



where x 12 = #i — x 2 , 



T-/ \ _ ( Q pi Q P2) 



and the quasiaverage 

/ rf 3 (x M (xi)rf 3 (T i ,(x2)p M p ly /(xi,p)/(x2,p)cos(ga;i2) 



(cos(ga?i2)> 



/ d 3 af J ,(x 1 )d 3 a u (x 2 )PiP2f(xi,pi)f(x2,P2) 
J d A X\d A X2 cos(qxu) S (xi,p) S (x 2 ,p) 



J d^xi S (xi,pi) J d 4 x 2 S(x2,p 2 ) ^ 
Obviously, the CF decreases from 2 at g = to unity at large |g M | and the width of this 
Bose-Einstein enhancement is inversely related to the effective system size which is usually 
encoded in so called interferometry radii. 

In this paper, we will also use Eqs. (JH)-© to calculate spectra and CF's at any earlier 
evolution times t a < t a on a hypersurface a. To do so, one has to make only the substi- 
tutions cr out — > a and S(x,p) — ► S a (x,p). Note that the current emission function S a (x,p) 
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takes into account the kinetic evolution up to the evolution time t a only and satisfies the 
obvious relations: S a (x,p) > S(x,p) for t < t a and S a (x,p) = for t > t a . Particularly, 
for a system of noninteracting particles the emission function is independent of the evolu- 
tion time and the asymptotic one, S(x,p), coincides with that, 5 ,<J °(a;,p), determined on the 
initial hypersurface a : t = to( x )- 

To calculate the kinetic evolution, we will use the simple ansatz S a °(x,p) = p f(x,p)5(t) 
for the initial emission function on the hypersurface o"o : t — 0. Though this ansatz 
generally violates the uncertainty relation At ■ Ap° > 1, it causes no problem since the 
evolution initial condition is determined by the initial distribution function only and since 
there exist a whole family of the initial emission functions of a finite time width recovering 
the latter in accordance with Eq. 

III. THEORETICAL ASPECTS OF THE SPECTRA FORMATION. 

The distribution function describes an expanding system of interacting particles and 
satisfies Boltzmann equation (BE) which, in the case of no external forces, has the form: 



The collision terms C ga i n and Ci OS s are associated with the numbers of particles that respec- 
tively come to the phase space point (,,,) and leave this point becanse of colhsions. Let 
us follow Ref. |l| and split the distribution function at each space-time point x into two 
terms: f(x,p) = f- m t(x,p) + f esc (x,p). The first one, fi Qt (x,p), describes the particles con- 
tinuing to interact after the time t. The second one, f eS c{x,p), describes the particles that 
do not interact any more after the time t. Introducing escape probability P(x,p), i.e., the 
probability that any given particle at phase space point (x,p) does not interact any more 
and propagates freely, one can write 



At large enough times t ~ t out , P(x,p) — > 1 and f eS ci x ,p) f( x ,p)- Since in the evolution 
of f esc (x,p) there is no loss of particles and the gain is associated with the part of the 
emission function, S g (x,p) = [1 — 9(t — t)]S(x,p), developed during the kinetic evolution 
starting at a hypersurface ctq and describing particles suffering last collisions at space-time 



P 



df(x,p) 
dx^ 



Cgain^j P) ~ C\ oss (x, p). 



(9) 




(10) 
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point x, the escape function satisfies the following BE 

d 

P^-r^fe S c{x,p) = P(x,p)C gain (x,p) = S g (x,p). (11) 

Noting that ^[p^exp (—iqx)] = and that f eS c(x,p) = f(x,p) on the asymptotic hypersur- 
face cr out , one obtains, applying the Gauss theorem to the Fourier transform of the distribu- 
tion function in Eq. (j3J) and using respectively general equations (|TT|) and analytically 
continued to off-mass-shell four-momenta p: 

r r a out 
«a P2 > = / d 3 a fl (x )p»f esc (x ,p)e-^+ / AS g (x )P )e^, (12) 



r r a out 
(a+a P2 )= / d 3 a ll (x )p»f(x ,p)e- iqxo + / d*x [C gaill (x, P )-C loaB (x,p)} e~ igx , (13) 

J (to «/crO 

where S g (x,p) = PC gain is defined in Eq. flTTj) . f eS c(%o,p) corresponds to the portion of 
the particles that are already free at initial time to (generally, at initial hypersurface cxq) 
and f(xo,p) is the distribution function at cr - Thus, as stated in Eqs. © and ©, the 
Fourier transform of the distribution or escape function on the asymptotic hypersurface cr out 
is equivalent to the four-volume Fourier transform of the emission function, the latter being 
represented by the emission function S g = PC gain developed during the kinetic evolution 
starting from the initial hypersurface o"o, together with the emission function 8(to — i)S(x,p) 
corresponding to emission times t < t UQ . The latter determines the escape function f esc (xo,p) 
on the hypersurface o"o through Eq. (JU) modified by the substitutions / — > f csc , S a — > 
6(t -i)S(x,p). 

One may see from Eq. (J 13)) that for a quasifree system characterized by the equality 
Cgain = Cioss, the spectra and correlations are not affected by the kinetic evolution and 
are completely determined by the initial distribution function. An example is the system of 
nonrelativistic particles with the initial spherically symmetric Gaussian distribution function 

DQ 

In hydrodynamic approach, the Landau freeze-out criterion of locally equilibrium (leq) 
hydrodynamic momentum spectra and corresponding Cooper-Frye prescription, defined by 
Eq. (jHJ) with substitutions cr ou t— > <J{ and /— > /i eq , treats particle spectra as a result of rapid 
conversion of a hadron system in local equilibrium into a gas of free particles on some freeze- 
out hypersurface Of. Formally, it corresponds to taking the cross section tending to infinity 



at t < t a{ (to keep the system in local equilibrium) and zero beyond t a{ . Then C ga i n = Ci OS s, 
P{x,p) = 6{t - t a{ ), / esc = 6{t - t CTf )/i eq and S = p°5(t - t CTf )/ieq- However, as stressed 
in Ref. [lo|, in the realistic case of continuous freeze-out the emission function S(x,p) is 
no more proportional to the distribution function f(x,p). At the same time, the spectra 
and the interferometry radii can be simply determined using the Cooper-Frye prescription 
for the fireball characterized by the thermal distribution function at the moment before it 
starts to emit particles and decays. As analyzed in detail in Ref. jl|, such a duality between 
hydrodynamic and transport approaches can happen when the integral of the difference 
Cgain — Cioss in Eq. (fT3|) vanishes, as is the case for a system of nonrelativistic particles with 
the initial spherically symmetric Gaussian distribution function (see sections IVI and rVT|) . or, 
when this integral is small, what can happen for a more wide class of initial conditions (see 
section IVII|) . 



IV. UNIVERSAL KINETIC MODEL. 

In our studies, we have used the UKM-R code Q written under ROOT system. In 
this code one deals with particle objects and their lists. Any particle object has its static 
and dynamic attributes. The static attributes are such as particle names, charges, masses 
or widths. The dynamic attributes are particle four-coordinates and four-momenta. The 
particle objects are grouped into particle lists. The ROOT classes are used to describe 
particle properties, to operate with three- and four-vectors and to organize containers. 

In the present work, we have considered only one type of particles and taken into account 
only their elastic binary collisions. We have thus created a specific list of primary particles 
specifying their spatial and momentum distributions and derived the Elastic interaction and 
Cascade model classes from the respective UKM Particle interaction and Universal kinetic 
algorithm ones. The Cascade model class represents the BE solver. The used UKM-R deals 
with the following lists: 

(i) Primaries, (ii) Collisions, (iii) Secondaries, (iv) New secondaries 
and the following basic classes: 

(i) Initial state, (ii) Elastic interaction, (iii) Cascade model. 

The numerical algorithm |(| can be described as follows: 
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1. The list of primaries is initialized and an initial value of the current time is set. 

2. At first, all possible pairs of colliding particles from the primaries list are searched 
for and the collision list is filled in according to the increasing collision time. The 
binary particle collision is considered as possible whenever the distance d of closest 
approach of the two particles in their cm. system is less than the collision distance 
calculated from the total collision cross section which is a function of the two-particle 
cm. energy y/s, i.e., 

d< [a(v/i)/7r] 1/2 . (14) 

In the present study, only elastic scattering is considered and a constant test cross 
section is used. 

3. The earliest collision is simulated, the two collided particles are removed from all 
lists and the created particles are put in the list of new secondaries. The particle 
coordinates are shifted in accordance with the time of simulated collision and particle 
velocity vectors. 

4. The collisions of the particles from new secondaries list with all particles from primaries 
and secondaries lists are searched for and the time ordered collision list is updated. 

5. The content of new secondaries list is moved to secondaries list and simulation contin- 
ues starting from step 3 or stops when the collision list is empty or the current time 
exceeds the user defined value of the stop time. 



It should be noted that the above algorithm violates the Lorentz covariance since 
the condition (fT4"|) introduces action at distance and breaks local character of the BE (see 
Eqs. (JOJ and ()15|1 ). One can use the invariance of the BE with respect to the transformation 
/ — > kf, a — > a/k and diminish nonlocality of the algorithm by simultaneously increasing 
the number of particles and the inverse interaction cross section by the same factor k. 

V. EXACT SOLUTION OF NONRELATIVISTIC BE. 

Let us consider the kinetic evolution of a system of N nonrelativistic identical spin-0 
particles of mass m assuming the initial conditions corresponding to the global thermal 
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Boltzmann distribution with temperature To (i.e., the absent initial flow) and spherically 
symmetric Gaussian profile of particle density with radius i?o- Assuming further only elastic 
binary collisions in the system, the collision integral in BE © takes the form 

C g ain(:r,p) -C loss (x,p) = J < ^^dn\f(x,p')f(x,p' 1 ) - f(x,p)f(x,pi))F(s)a(y/s,6), (15) 

where the particle four-momenta in the binary scattering satisfy the conservation law: p + 
Pi = p' + Pi, dVt = dcos8d<f), cos# = 1 + 2t/(s — 4m 2 ), s = (p + P\) 2 , t = (p — p') 2 , 
F(s) = — 4m 2 )] 1//2 is the invariant flow of two colliding particles of mass m at cm. 

energy ^/s and <j(^/~s,9) = da e i/dQ is the differential cross section of the elastic scattering 
on the cm. angle 9. Then, in the nonrelativistic case, the BE has the following analytical 
solution for the distribution function f(x,p) = f(t, r, p) Q]: 

fit \ N ( (r - tp/mf p 2 \ 
/(i ' r ' P)= |2^(»T )Vf eXP ( M SSi.l- <16) 

Since Eq. f!16|) represents also special case of the solution of hydrodynamic equations for 
ideal fluid with the radial initial flow, it can be rewritten in the local equilibrium form jjj: 

N f r 2 (p-mu(t,r)) 2 ' 

At, r, p) - /lcq - [27ri?o(mTo)1/2]3 exp ^ 2i?2(t) 2mT{t) 

where T(t) = T /(p(t), R 2 {t) = B%(p(t) and u(t, r) = T(p(t)/(p(t) = vtT(t)/(mR 2 ) are 
respectively the temperature, the system Gaussian radius squared and the flow velocity, 
ip{t) = 1 + T t 2 /(mRl); the corresponding chemical potential (i(t,r) = —^T{t)r 2 /R 2 {t). 

Inserting solution (fT7)|) into the BE, one can check that its l.h.s. vanishes, mdf/dt + 
pdf/dr = 0, so the collision integral (its r.h.s.) should also vanish. Indeed, inserting the 
solution (J 17)) into the factor in square brackets in Eq. (|15|) . one can see that the collision 
integral vanishes (irrespective of energy and angular dependence of the elastic cross section) 
as a consequence of the energy and three-momentum conservation: 

[f(x,p')f(x,p[) - f(x,p)f(x,p 1 )] oc exp (- P tP* + P -^ u(t,r) 



2mT(t) T(t) 
( P 2 + Pi , P + Pi ( . C\ n na \ 

This solution thus corresponds to a quasifree motion and so, as discussed in section EH (see, 



e.g., Eq. (O), the particle spectra and 
and stay unchanged at any times 



correlations are not affected by the kinetic evolution 
One may see from Eqs. (fTBj) and (|T7)) that the 



conservation of the momentum spectrum is related to the fact that the decreasing temper- 
ature is compensated by increasing flow velocity. The latter also compensates the effect of 
increasing system size on the two-boson CF thus guaranteeing the conserved interferometry 
radius. Indeed, using distribution function explicit calculations yield time-independent 
momentum spectrum and correlation function: 



d 3 N(t; p) 
d 3 p 



N 



(2vrmT ) 3/2 



exp 



P^ 



2mTn 



(19) 



ft(*;Pi,P2) = 1 + exp 



R l 



4mTn 



(20) 



where (4mT )^ 1/2 = R th measures the characteristic thermal size (heat de Broglie length) of 
the single-pion emitter. At the same time, as seen from Eq. (|17j) . the system expands and 
spatial density decreases with increasing time: 



d 3 N(t; r 



N 



d 3 r {2^{t)f 2 ^ V V 1R\t)) 



(21) 



where R 2 (t) = R 2 + t 2 T /m. 



VI. TEST OF THE UKM-R ALGORITHM. 



To test the UKM-R algorithm, we have implemented the above spherically symmetric 
Gaussian initial conditions and simulated the kinetic evolution in a number of events with 
fixed boson multiplicity N = 400, choosing To = 0.130 GeV and Rq = 7 fm. To guarantee the 
nonrelativistic particle momenta, we have put m = 0.938 GeV/c 2 . The initial momentum 
and coordinate distributions are shown in Fig. ^ We have done two sets of simulations 
with the elastic cross sections a e i = 40 and 400 mb, assuming the isotropic scattering in the 
two-particle cm. system. We have checked that the final results remain unchanged when 
introducing the cross section anisotropy. 

Instead of the six- dimensional CF in Eq. ©, we have calculated one-dimensional ones 
integrating in numerator and denominator over all kinematic variables except a chosen one, 
e.g., the modulus q = |q| of the relative momentum vector or its projections, or, the invariant 
relative momentum q inv = (q 2 — Qq) 1 ^ 2 = q(l — v 2 ) 1 / 2 , where v g is the projection of the pair 
velocity vector on the q-direction; for nonrelativistic particles ^ C 1 and qi nv = q. As usual, 
we have assumed sufficiently smooth behaviour of the single-particle momentum spectrum 
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(i.e., in our case, the initial thermal radius R t h much smaller than the initial system radius 
Rq) and substituted the quasiaverages in Eq. (jSJ) by the averages making the respective 
substitutions x iy p — > x i: pi and Xi,p — > Xi,Pi in the arguments of the distribution and 
emission functions in the numerators. Explicit calculation at t — then leads to the Gaussian 
CF in Eq. (J2U|) with the substitution R^ — R% b — > R^. It should be noted however that this 
approximation violates the invariance of the four-vector product qx on the trajectory of a 
free streaming particle (the invariance is guaranteed only for free streeming particle with the 
off-mass-shell four-momentum p = (pi +p 2 )/2). As a result, the averages with distribution 
and emission functions are no more equivalent, the latter being preferable since it avoids 
free streaming from the emission or collision point to a given hypersurface a. 

Concerning the approximate averaging using the distribution function, at large evolution 
times t as compared with the the mean emission time (t), it leads to the additional oscillating 
factor in the CF, 

TZ{t; Pl , p 2 ) = 1 + exp (-i? V) cos (J^qLj ■ (22) 

The oscillating factor invalidates the smoothness approximation at very large times and, at 
moderate times, it leads to an increase of the correlation radius, roughly, 

%-)', (23) 

where the factor k ~ 1. For heavy particles of mass m = 0.938 GeV/c 2 and the initial 
system radius R = 7 fm, this increase does not exceed 10% provided t < 100 fm/c. 

Similarly, the approximate averaging using the emission function at a current evolution 
time t yields 

TZ(t; Pl , p 2 ) = 1 + exp (-R-W) (cos (—q? nv ) ) , (24) 

\ \ m J I tl, S <t 

where i\2 = (ti + ^)/2 and the averaging is done over the emission times ti^ < t. For 
large evolution times, t ~ t out , the distortion factor in Eq. (|2^|) is independent of t and its 
effect on the CF is determined by (t^) t out . It thus leads to only a slight increase of the 
correlation radius, roughly determined by Eq. (|2l?j) with the substitution t — > (i^). 

Practically, we have calculated the one-dimensional CF in the bins of a chosen variable 
q j (e.g., q inv ) as 

AN(q 3 ) 
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where AN(qj) is the number of simulated pairs in a given bin, consisting of the particles 
from the same event only; i.e., each pair from the same event was attributed the weight [1 + 
cos(q^Xil)] and the weighted histogram of qj was produced. To determine the corresponding 
interferometry radius Rj, the constructed CF has been fitted as 

K{q j ) = l + \exp(-B$<ij)- (26) 

To check the exact solutions for the momentum and coordinate distributions in Eqs. (jl9j) 
and (j21j) and the CF in Eq. (|2"U|) . the particle coordinates and momenta have been evaluated 
at the following values of the evolution time t: 0, 50 and 100 fm/c. In Fig. [I]we show the 
results of the UKM-R simulations corresponding to the elastic cross section o e \ = 400 mb; we 
have checked that they coincide with those obtained for a e \ = 40 mb. These results confirm 
the independence of the momentum spectrum and the CF on the evolution time despite that 
a huge number of collisions happened during the evolution. A slight increase of the fitted 
interferometry radius at t = 100 fm/c is a consequence of the smoothness approximation 
and agrees with Eqs. ()22|) and ()23|) . As for the coordinate distribution, panel b of Fig. ^ 
confirms the increase of the system Gaussian radius with the evolution time in accordance 
with the law R 2 {t) = R% + t 2 To/m. The results shown in Fig. Q]are in full correspondence 
with the exact solution of the BE for distribution function and thus confirm the correctness 
of the kinetic simulation. 

As noted after Eq. (}2"2"|) . the use of the smoothness approximation in the calculation of 
the CF with the help of the distribution function is invalidated at large evolution times. 
This is demonstrated in Fig. |21 where the appearance of the artificial zero of TZ(qi nv ) — 1 
at qi nv = (irm/2t) 1 /' 2 « 0.02 GeV/c due to the oscillating factor in Eq. (|22|) is clearly 
seen in the CF calculated from boson coordinates at the evolution time t = 700 fm/c. 
To calculate the correct CF at asymptotic times, we have therefore followed the system 
evolution up to a sufficiently large time (t ~ 200 fm/c for a e i = 400 mb) and taken as 
the particle four-coordinates those of the last emission or collision points. This corresponds 
to the system description with the help of the asymptotic emission function S(x,p) which 
is much less affected by the smoothness approximation as compared with the distribution 
function. The spatial and time distributions of these points, corresponding to a e \ = 40 and 
400 mb, are shown in the panels b and c of Fig. El For a e i = 40 mb, one may see only a 
slight broadening of the spatial distribution as compared with the initial one and a quite 
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narrow time distribution with a noticeable number of nonscattered particles (t = 0) and 
the asymptotic time of ~ 100 fm/c. For o e \ = 400 mb, the spatial and time distributions 
are much broader, almost all particles scatter at least once and the asymptotic time is 
~ 200 fm/c. As for the asymptotic momentum spectra and interferometry radii, one may 
see from the panels a and d of Fig. |3] that they coincide with the initial ones within the 
errors for both values of a e i. 

VII. ANISOTROPIC, NON-GAUSSIAN AND RELATIVISTIC SYSTEMS. 

We have also studied to what extent are the spectra and CF's conserved in kinetic evo- 
lution of nonrelativistic and relativistic systems characterized by initial conditions different 
from the spherically symmetric Gaussian ones given in Eq. (JTfij) at t — 0. 

A. Nonrelativistic systems. 

1. Anisotropic non-Gaussian spatial distribution. 

As in section IVI| we have considered the nonrelativistic system of N = 400 heavy spin-0 
bosons of mass m = 0.938 GeV/c 2 and the initial momentum distribution given in Eq. (fTTjj) 
with T = 0.130 GeV. However, the initial spherically symmetric Gaussian spatial distri- 
bution was substituted by a uniformly populated parallelepiped defined by \x\ < 14 fm 
and \y\, \z\ < 7 fm (Fig. 0j); the corresponding initial volume is about twice as large than 
in the previous case. The kinetic evolution was performed with two values of the elastic 
cross section, a e \ = 40 and 400 mb. In the former case, the collisions are quite rare so the 
original spatial and momentum distributions, as well as the CF, are practically not changed 
during the evolution. In the latter case, the collisions are much more frequent and force the 
spatial distribution of the emission or last collision points to become close to a spherically 
symmetric Gaussian one (see dashed histograms in panels a and c of Fig. HJ). On the other 
hand, as shown in Fig. |5J the initial space anisotropy is transformed in the anisotropy of the 
final three-momentum distribution. As compared with the originally spherically symmetric 
three-momentum distribution, the final one becomes softer in the x-direction and harder 
in the (y, z)-plane. As for the CF's, one may see from Fig. |U] that their evolution reflects 
the one of the emission points, the final CF's becoming closer to spherical symmetry. It is 
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interesting to note that despite of the noticeable difference of the initial and final CF's, the 
interferometry volume, R X R^, is practically not changed during the evolution, the initial and 
final one composing 177 ±6 and 182 ±6 fm 3 respectively. Also the qi nv projection of the CF 
is changed only slightly (see panel c in Fig. The effect of approximate conservation of 
the interferometry volume is in correspondence with a similar result found for the chemically 
frozen hydrodynamic evolution 



2. Nonthermal momentum distribution 

We have further considered the same initial conditions as in section IVII except for the 
thermal Gaussian three-momentum distribution which has been substituted by a uniformly 
populated cube defined by \p x \, \p y \, \p z \ < (T m) 1//2 , T = 0.130 GeV, m =0.938 GeV/c 2 
(Fig. |7J). The results of the kinetic evolution of the momentum distribution and the CF's 
corresponding to the initial Gaussian radius Ro = 7 fm and the elastic cross section a e \ = 400 
mb are shown in figures and El similar results have been obtained with a e \ = 40 mb. One 
may see that the collisions force the initially nonthermal distribution to take the thermal 
Gaussian form. At the same time, r.m.s. of the initial distribution of the momentum 
components, a = (mT /3) 1//2 = 201.6 MeV/c, as well as the initial CF and the interferometry 
radius of 7 fm, are practically conserved. E.g., the Gaussian fits of the p x distribution and 
the q t projection of the CF yield a = 202.4±0.2 MeV/c and R t = 7.14±0.11 fm respectively. 



B. Relativistic systems. 

We have also studied the kinetic evolution of the relativistic hadronic gas (m ~ T ) 
choosing the particle mass m = 0.140 GeV/c 2 and the initial temperature To = 0.130 GeV. 
It is well known that the only relativistic solution of the BE with vanishing collision term is 
the ^-independent Jiittner distribution function, 

/M=(4exp(^), (27) 

corresponding to a global thermal equilibrium with temperature T, chemical potential fi 
and system four-velocity U . Particularly, the relativistic generalization of the distribution 
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function in Eq. (fTol) . 

f{x,p)ocexp^ —J, (28) 

represents only an approximate solution of the BE, leading to vanishing l.h.s. but retain- 
ing a nonzero r.h.s. collision term. Generally, one may expect that, irrespective of initial 
conditions, the distribution function formed in the kinetic evolution should tend to the 
Jiittner one, i.e., the spatial distribution should expand and the three-momentum distribu- 
tion d 3 N/d 3 p should tend to the thermal one cx exp[— (p 2 + m 2 ) 1 / 2 /T]. Particularly, starting 
the kinetic evolution from the initial conditions corresponding to Eq. ()28j) at t — 0, one may 
expect an increase of the system Gaussian radius with the evolution time and conserved 
momentum distribution. Based on the analogy with the nonrelativistic system, one can also 
expect an approximate conservation of the correlation radius. To check these expectations, 
we have put the initial Gaussian radius Ro = 7 fm and simulated the kinetic evolution 
for N = 400 particles with the elastic cross section cr e / = 400 mb. The simulation results 
shown in figures E] and El indeed confirm the conservation of the momentum spectrum and 
approximate conservation of the CF, the final correlation radius exceeding the initial one 
by ~ 9% only. We have obtained the same results for N = 40000 particles and a e \ = 4 mb 
thus proving the negligible effect of the violation of Lorentz covariance introduced by the 
nonlocal condition in Eq. (fPfj) ; some deviations start to appear only at a e i ~ 1000 mb. 



1. Anisotropic non-Gaussian spatial distribution. 

We have considered the same initial anisotropic coordinate distributions as for the nonrel- 
ativistic bosons. Comparing figures HHBI and ITTll 1 31 corresponding respectively to nonrelativis- 
tic and relativistic bosons, one observes similar isotropisation effects in the distributions of 
the emission points and CF's, as well as similar anisotropy transfer from spatial to momen- 
tum components of the phase space emission points. The interferometry volume is however 
conserved only approximately, it increases by ~ 30% from the initial value of ~ 176 fm 3 to 
the final one of ~ 250 fm 3 . 
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2. Nonthermal anisotropic momentum distribution 

We have further considered the same initial conditions as in Fig. |§] except for the thermal 
momentum distribution which has been substituted by a uniformly populated parallelepiped 
defined by \p x \ < 0.7 GeV/c, \p y \, \p z \ < 0.35 GeV (Fig. Hljl . The results of the kinetic 
evolution of the three-momentum distribution and the CF's are shown in figures ITU and [TBI 
One may see that the collisions force the initially nonthermal distribution to become close 
to the thermal one. The final interferometry radii are ~ 10% higher than the initial one of 
7 fm, similar to the case of the initial thermal momentum distribution. 



VIII. CONCLUSIONS AND OUTLOOK 



We have studied numerical solutions of the Boltzmann equation in the Universal Kinetic 
Model. First, we have tested the model kinetic algorithm recovering the known nonrel- 
ativistic solution [1-3] for the spherically symmetric Gaussian initial conditions up to the 
elastic cross section of 1000 mb. Particularly, we have checked that despite a huge number of 
collisions the momentum spectra and interferometry radii remain unchanged during the evo- 
lution. Further, we have modified the initial conditions and found that with the increasing 
elastic cross section the system of nonrelativistic particles more and more recovers spherical 
symmetry and thermal momentum distribution. At the same time, similar to the above 
case, the momentum dispersion and the interferometry volume is practically conserved dur- 
ing the evolution. The conservation effect takes place also in the evolution of the system of 
relativistic particles, except for some increase (~ 30% at 400 mb) of the final interferometry 
volume. Our studies thus support similar results obtained in hydrodynamic approximation 
inRef. Q. 

It is worth to note that the initially anisotropic system tends to an isotropic one and, 
after sufficiently large number of collisions, the spectra and interferometry radii are approxi- 
mately conserved in subsequent evolution. In such a case one can exploit the idea of duality 
of hydrodynamic and kinetic de scrip tions of interferometry volumes and effective temper- 
atures (momentum dispersions) [lOj, i.e. estimate their values from the thermal hadronic 
distribution function in the initially formed hadronic fireball instead of a detailed analysis 
based on the complicated asymptotic distribution or emission functions. 
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We plan to extend the analysis of kinetic evolution of the spectra and interferometry radii 
to many component hadron gases, including inelastic collisions and resonance decays. 
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FIG. 1: Distributions of the x-components of particle three- momentum (a) and three-coordinate 
(b), particle momentum (c) and the corresponding CF's as functions of qi nv (d) obtained in the 
UKM-R simulation of the kinetic evolution of N = 400 heavy spin-0 bosons of mass m = 0.938 
GeV/c 2 at the evolution time t = 0, 50 and 100 fm/c. The elastic cross section a e i = 400 mb, the 
initial Gaussian radius Rq = 7 fm and the initial temperature To = 0.130 GeV. The shown results 
of Gaussian fits of the momentum distributions and the CF's agree with the input initial values 
a = (mTo) 1 / 2 = 0.3492 GeV/c and A = 1, Rq = 7 fm respectively. 
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FIG. 2: CF's obtained from four-coordinates of free streaming {a e i = 0) heavy spin-0 bosons of 
mass m = 0.938 GeV/c 2 obtained in the UKM-R simulation on the same conditions as in Fig. Q 
at the evolution times 100 and 700 fm/c. 
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FIG. 3: Initial and final distributions of the x-components of particle three-momentum (a), spatial 
x- (b) and time (c) coordinates of the emission or last collision points and the CF's corresponding 
to these points as functions of qi nv (d) obtained in the UKM-R simulation of N = 400 heavy spin-0 
bosons of mass m = 0.938 GeV/c 2 at the evolution time t = 200 fm/c. The histo grams 1, 2 and 3 
in panels b and c correspond to a e i = (initial x-distribution), 40 and 400 mb respectively. The 
shown results of Gaussian fits of the momentum distributions and the CF's agree with the input 
initial values a = (mTo) 1 / 2 = 0.3492 GeV/c and A = 1, R = 7 fm respectively. 
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FIG. 4: Final distributions of x- and y-coordinates of the emission or last collision points (dashed 
histograms) obtained in the UKM-R simulation on the same conditions as in Fig. |3l except for the 
initial distributions of spatial coordinates (solid histograms): |x| < 14 Fm, |y|,|z| < 7 Fm. The 
panels a, c and b, d correspond to a e \ = 400 and 40 mb respectively. 



22 




Px, GeV/c 




FIG. 5: Initial (thick) and final (thin) histograms of the components of particle three-momentum 
obtained in the UKM-R simulation on the same conditions as in Fig. 21 
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FIG. 6: Initial and final CF's corresponding to the UKM-R coordinate distributions in panels a 
and c of Fig. |I] (a e i = 400 mb). The CF's are shown as functions of q x (qt < 10 MeV/c) (a), 
Qt = (Qy 2 + Qz 2 ) 1 ^ 2 (Qx < 7 MeV/c) (b) and qi nv (c). The results of Gaussian fits are shown. The 
fitted initial Gaussian radii in panels a and b are close to the r.m.s. values 8.08 and 4.04 fm of the 
corresponding initial x- and y-coordinate distributions in panels a and c of Fig. 21 
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FIG. 7: Final distributions of the x-component (b) and the magnitude of particle three-momentum 
(d) obtained in the UKM-R simulation with a e i = 400 mb on the same conditions as in Fig. |3] 
except for the initial three-momentum distribution (a and c) defined by \p x \, \Py\, \Pz\ < (Tom) 1 / 2 , 
Tq = 0.130 GeV, m =0.938 GeV/c 2 . The fitted Gaussian widths of final distributions, shown in 
panels b and d, are close to the r.m.s. value (mTo/3) 1 / 2 = 0.2016 GeV of the initial p x -distribution. 
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FIG. 8: Initial and final CF's corresponding to the UKM-R simulation specified in the caption 
of Fig. El The CF's are shown as functions of q x (q t < 10 MeV/c) (a), q t = (q y 2 + q z 2 ) l/2 
(q x < 7 MeV/c) (b) and qi nv (c). The results of Gaussian fits are shown. The fitted correlation 
radii are close to the initial system Gaussian radius Rq = 7 fm. 
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FIG. 9: Distributions of particle energy (a, c), momentum (d) and x-component of particle 
three-momentum (b) obtained in the UKM-R simulation of the kinetic evolution of N = 400 light 
spin-0 bosons of mass m = 0.140 GeV/c 2 and elastic cross section 400 mb. The initial three- 
coordinates and three-momenta have been respectively distributed according to Gaussian of the 
radius Rq = 7 fm and relativistic thermal distribution with the temperature To = 0.130 GeV. The 
initial distributions are shown as solid histograms in panels a, b, d and the final ones - as solid 
histogram in panel c and dashed histograms in panels b, d. The fitted values of the initial and 
final temperature shown in panels a and c coincide with Tq. 
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FIG. 10: Initial and final CF's corresponding to the UKM-R simulation specified in the caption 
of Fig. The CF's are shown as functions of qi nv (a) and q = |q| (b). 
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FIG. 11: Final distributions of x- (a) and y- (b) coordinates of the emission or last collision points 
(dashed histograms) obtained in the UKM-R simulation on the same conditions as in Fig. |2] except 
for the initial distributions of spatial coordinates (solid histograms): \x\ < 14 Fm, |y|, \z\ < 7 Fm. 
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FIG. 12: Initial (solid) and final (dashed) histograms of the components of particle three- 
momentum obtained in the UKM-R simulation on the same conditions as in Fig. 1111 
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FIG. 13: Initial and final CF's corresponding to the UKM-R simulation specified in the caption of 
Fig.HTJ The CF's are shown as functions of q x (q t < 10 MeV/c), q t = (q y 2 + q z 2 ) 1/2 (q x < 7 MeV/c) 
and q = |q|. 
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FIG. 14: Final distributions of the components (b, d) and the magnitude (f) of particle three- 
momentum obtained in the UKM-R simulation with a e i = 400 mb on the same conditions as in 
Fig. 03 except for the initial three-momentum distribution (a, c, e) defined by \p x \ < 0.7 GeV/c, 
\Py\i \Pz\ < 0.35 GeV/c. The fitted Gaussian widths of final distributions are shown in panels b, d 
and f. 
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FIG. 15: Initial and final CF's corresponding to the UKM-R simulation specified in the caption of 
Fig.Qll The CF's are shown as functions of q x (q t < 10 MeV/c), q t = (q y 2 + q z 2 ) 1/2 (q x < 7 MeV/c) 
and q = |q|. 
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